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© Frequency control apparatus and method for a digital radio receiver. 

© A frequency control apparatus for a burst-mode radio communications system (TDMA) employina a 
ET* C f or : eCt,0 :. S,9na, ^ a burst to enable correction of frequency differences bXeen the 

23*t£ * C3mer ° f 3 SUbSequent communication burst and the frequency of iSto^JE 

receiver, the frequency control apparatus including a variable frequency oscillator (115) fat receivino a 
frequency correction signal burst and a digital signal processor (111) for processing I c^iS^ri 
having a frequency difference between said variable frequency oscillator and the radio c^ 
communication burst and for generating a singular control signal value and for eJy^Zid^Z clnt^ 
s.gna, value to said variable frequency oscillator (115). via an converter, to corSS 
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FREQUENCY CONTROL APPARATUS AND METHOD FOR A DIGITAL RADIO RECEIVER 



Background of the Invention 

This invention relates generally to digital radio systems and more particularly to the rapid comcton of 
frequency error in a digital radio rece.ver which rece.ves communication messages from a transm.tter 
s transmitting the communication messages in bursts. 

Radio receivers often utilize a c.rcuit to automatically correct for discrepances ,n frequency between 
the carrier frequency of the s.gnal to be received and the frequency of a local oscllat or used m a 
superheterodyne receiver. This local osc.llator converts the earner frequency and associated >" °™*mn 
carrying sidebands of the received signal to a convenient intermediate frequency. The typical frequency 
,o correction process is performed over a relatively long period of time, assuming that the earner ***** £ 
the received signal) is continuously present. The carrier frequency may be tracked by means of amplitude 
detectors, discriminators, or the like to generate a frequency control signal. Some systems may use a pilot 
s.gnal modulated onto the earner frequency to provide a reference in the receiver to derive the frequency 
control signal. (Such a pilot control is further described in U.S. Patent No. 4.541 .118). 
, 5 The frequency control signal is subsequently applied to the local oscillator from its initial frequency to a 
frequency which converts the received earner frequency into an intermediate frequency optimally placed 
within the selectivity of the intermediate frequency amplification and filtering stages. 

Digital receivers present a new set of problems to such conventional automatic frequency control 
networks. One approach to a fast receiver frequency control for a digital receiver has been d.sc osed ,n 
20 European Patent Application No. 89122616.9. "Digital Automatic Frequency Control on Pure Sine Waves . 
filed on 11th December 1988 on behalf of Motorola Inc. 

Usually, digital receivers must process the received carrier frequency signal in a linear fash.on. Such 
linear processing allows amplitude variation of the received signal to create further errors in the detection of 
frequency offset. Furthermore, digital communication is often accomplished using burst transmission 
" 25 techniques such as time division multiple access (TOMA). Burst transmissions do not prov.de a continu- 
ously available carrier or carrier plus pilot which can be employed by conventional frequency control 
techniques. 



30 Summary of the Invention 

Therefore, it is one object of the present invention to rapidly correct for frequency errors between a 
received signal and the frequency to which the receiver is tuned in a digital radio receiver. 

It is another object of the present invention to rapidly correct a receiver oscillator for frequency errors 
35 between a nominal frequency and a desired frequency. 

It is another object of the present .nvention is to correct the frequency error between the received signal 
and the correct receiver local oscillator frequency within a single frame of received data. 

It is a further object of the present invention to determine the frequency error from a correction burst 
and directly correct the local oscillator frequency. 
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Brief Description of the Drawings 

Figure 1 is a block diagram of a TDMA receiver which may utilize the present invention. 
Figure 2 shows the unwrapping process for both an increasing phase trajectory and a decreasing 
phase trajectory. 

Figure 3 shows the result of a computer simulation. 

Figure 4 shows a timing diagram relating the burst TDMA transmissions to the frequency control of 
the present invention. 

Detailed Description of the Preferred Embodiment 

A block diagram of a TDMA receiver which may utilize the present invention is shown in Fig. 1. Such a 
receiver may receive burst communication signals in a timeslot from a TDMA transm.tter 103 received on 
an antenna 104 and applied to quadrature demodulator 105 as a signal x(t). The quadrature demodulator 
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105 produces two quadraturely related downconverted signals which are applied to conventional analog to 
digital converters 107 and 109 which digitize each of the two quadrature related signals before applying the 
digitized quadrature signals to a digital signal processor (DSP) function 111. In the preferred embodiment 
the DSP is realized employing a DSP56001 available from Motorola, Inc. (or equivalent). The DSP function 
1 1 1 recovers a data communication signal using conventional TDMA techniques in a communication signal 
recovery function 113. Such signal recovery functions include channel estimation, channel equalization and 
data demodulation. Also included in the DSP function 111 is the recovery of a frequency control signal 
which allows the rap«d adjustment of a voltage controlled oscillator 115 in the quadrature demodulator 105 
so that the TDMA receiver can be quickly brought into a condition where there is virtually no frequency 
error between the carrier frequency of x(t) and the frequency of voltage controlled oscillator 115. Because 
speed of adjustment is important . it is an important feature of the present invention that the adjustment of 
the frequency of voltage controlled oscillator 115 be substantially accomplished in one step without 
oscillator frequency hunting or successive approximations to a final oscillator frequency. 

,r,o,?,f Ume th3t 018 transmitted s, '9na' corresponding to a frequency correction burst (or pure sine wave 
(PSW)) is given by 

x(t) = A cos[« c t + « 0 t] (1) 

where: 

A is the amplitude of the signal. 

«c = 2irf c is the carrier frequency of the signal in radians, sec: 
o> 0 = 2wf 0 is the PSW baseband signal frequency in radians/sec 
= 2* (67.708...kHz). 

In the absence of multipath noise. Rayleigh fading, etc.. the received signal is given by (1). In general 
the voltage controlled oscillator 115 (VCXO) frequency without frequency control will be given by («* + A J 
where &<* corresponds to a frequency offset error (in radians/sec). AssurnTthTt the VCXO output signal is 
therefore given by 
v,(t) a cos[(« c + A«)t]. (2) 

After phase shifting v,(t) by -90* in the phase shifter 117. the resulting signal is given by 
v Q (t) = sin[(« c + Aa,)t]. (3) 

The output of the in-phase (I) mixer 1 1 7 (in the absence of any DC offsets) is given by 
IaW - x(t) v,(t) = A cos [« c t + « 0 t] cos[(o, c + A u )t] = 1/2 A { cos{(« 0 - A«)t] + cos [(2« e t + « 0 + A<-)t]} 

rH!° h " f o r J 0VV ^f SS fi ' terin9 With fHter 121 having a bandwi dth of approximately 4 radians/sec. becomes 
l B (t) = 1/2 A cos[(u 0 - A«)t]. (5) 

^ Similarly the output of the quadrature-phase (Q) mixer 123 (in the absence of any DC offsets) is given 
Qe(t) = x(t) v Q (t) = A cos [« c t + Mo t] S in[(*, e + A «)t] = 1/2 A {sin [(2o, c t + + A*)t] - sin [(«„ - A«)t]} 

which after low-pass filtering by filter 125 becomes 
Q B (t) = -1/2 A Sin[(<j 0 - Au)t]. (7) , 

The two A/D converters (107. 109) convert l 8 (t) and Q e (t) into their quantized, sampled equivalents. I 8 (k) 
and Q 8 (k). respectively. In the absence of any DC offsets: l B (k) = l c (k) and Q a (k) = Q c (k) Note that direct 
measurement of the frequency offset A« via a frequency counter which counts zero crossings of the signals 
given in (5) and (7) will yield a measurement uncertainty of 
freq. uncertainty = 1 /measurement period = 1/slot period 
= 1/0.58 msec = t 1724 Hz (8) 

which corresponds to an uncertainty of * 1.9 Ppm at 900 MHz. In certain radio systems, for example, a 
digital radiotelephone system for use in Europe defined in GSM documents GSM 05.01/3.1 0 "Physical 
Layer on the Radio Path: General Description" and GSM 5.10/3. "Radio Sub-System Synchronization", the 
required stab, ity of the VCXO 115 must be within * 0.1 ppm of the received signal coming from the base 
station^ Therefore, an approach to measuring frequency errors which does not count zero-crossings must 
instead be employed. The apparatus and method of the present invention- 

presenUn^a'nd 6 ^)? 6 ' ° S ' 9na ' S 9 ' Ven ^ U) ^ ^ kH2 l ° ° C * * ' rem ° V8S the "° term 

(2) computes the phase of the translated signal at regularly spaced Hme intervals: 

(3) constructs the phase trajectory from the phase samples: and 

(4) computes an estimate of the instantaneous frequency of the frequency-translated signals from the 
time-derivative of the phase trajectory. 
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DC Offset Compensation 

The DSP function 111 accepts the quantized l*Q signals and processes both the I and the Q signal 
s through a direct current (DC) subtraction process 129. An uncorrected DC offset is a result of mismatched I 
and Q channels and local oscillator leakage. An estimate of the DC offset of l s (k) is made by DC estimate 
function 131 which averages the input signal the output of which is subtracted from i B (k) in digital adder 
133 thereby producing l c (k). the I quadrature signal with the DC removed. Likewise, an estimate of the 
direct current offset is made by DC estimate function 135 of the Q B fk) signal and subtracted in digital 
jo adder 137 to produce Q c (k). the Q quadrature signal with the DC removed. 

In order to implement the above frequency control of the present invention in a fractional fixed point 
general purpose digital signal processor (DSP) such as the DSP56001. it is necessary to perform several 
unique steps. Imperfections in an implementation of the quadrature demodulator 105 necessarily results in 
DC offsets at the output of the quadrature demodulator 105 which must be compensated for in the 
rs frequency control. 

The presence of DC offsets at the output of the quadrature demodulator 105 can severely limit the 

performance of the frequency control. To see why this is the case, let l B (t) and Q 8 (t) in (5) and (7) contain 

DC offset terms of magnitude C and D, respectively. Then 

l Q '(k) = 1 2 A cos [<a> 0 - A«)k] + C (5a) 
20 Q B fk) a -1 2 A sin [(w 0 - Aa»)k] «• D (7a) 

Then after frequency translation by e ,u> ° k as described below 

l 0 "(k) = Re {[S B fk) + C + jD] eJ w ° k } 

= Re {[1 2 A e'* w ° ' Aw)k + C + jD] e* w ° k } 

= 1 2 A cos Au)k + C cos <j 0 k - D sin « 
25 Qo'lk) = im £[S 8 (k> + C + jD] e>» ° k } 

= lm{[1 2Ae"° , ° Au,)k + C + jDJe^Q *} 

= 1 2 A sin Awk + C sin *> 0 k + D cos w 0 k. 

The presence of the two quadrature components at a frequency w 0 prevent further downsampling (or 

decimation) without introducing aliasing errors and can contribute to a significant error in the computation of 
30 tan* 1 (Q E (k) le(k> ) in (13) below. Removal of the two undesired frequency components in (17) below can 

be accomplished by two methods: 

(1) Low-pass filtering of l D 0o and Qo'fk) to remove the sin « 0 k and cos o> Q k components; or 

(2) Removal of the DC components prior to frequency translation by e iQ> ° k . 

The first approach requires a low-pass filter with a bandwidth less than f 0 /10= 6.7kHz. for example, and 
35 thus having an impulse response of duration greater than 5/6.7 kHz = 0.73 msec, for example, an impulse 
response duration greater than the duration of a single time slot. Clearly this is an unacceptable approach. 

The second approach simply computes the average DC value of the two quadrature branches 
separately and subtracts these values from the respective branch signals: 

J 

I C (k) = I B (k) - 1/J X lB<» 

k=l 
J 

QC<k) = QbOO - 1/J X QB(k) 

k=l 



Frequency translation by 67.7 KHz. 



A TDMA receiver compatible with GSM specifications such as those mentioned above utilizes a 
frequency correction signal which is transmitted as a carrier shift of exactly 67.7 kHz during one timeslot of 
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the TDMA transmission. A frequency translation of the I quadrature channel signal and the Q channel 
quadrature signal by 677 kHz accomplishes the following result. The signal is translated to 0 Hz * Z7l 
thereby reducing the information data rate and permitting further decimination 
s 1 M ^ * a f omp,ished ^ a quantized 67.7 KHz signal in quadrature, as generated by 67.7 KHz oscillator 
r!L^TZ°OSP. aCCOmP " Shed by COmp,e * ^ adratu ' e ™>™ ^ich is a complex multiplication 
Let 

S c (k> = l c (k) + jQ c (k) = 1 2 Ae «"o - Aa>)k (9) 

where l c (k) and Q c (k) are given above. After multiplication of S c (k) by e» u ° k we obtain 
>o So(k) 3 , 0(k ) + jQo (k, = 1.2 A{e^*} - 1.2 A cosA.k + j 12 A sinLk '(10) 
i.e.. 

lo(k) = 1/2 A cos A«k 

Qo(k) = 12AsinA«k. (11) 

In the preferred embodiment, the TDMA modulation is conventionally processed by the transmitter in a 
' 5 oenTrSv r n "J ^ ?** ^ ° CCUPied bandwidth ° f the Emitted sTncT 2 win 

Z a norma "u sH^ead oTff inf ° rmati ° n ° f ^ GMSK Si 9" a * with -iom dlta 

*1 cl T frequency correction burst). S D (k) may be decimated by a factor of M to 

reduce the s.gnal processing complexity required for frequency control without reducing contro accuracy 
The decmated signal S E (k) is obtained from S 0 (k) by the operation accuracy. 
20 S E (k) = So(Mk) (12) 
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Computation of the phase of the decimated, translated signal. 
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vrvn TOT?*, tri S° nom8tric ^entities, a calculation of the phase relationship between the 

S ' 9n ^ 3 u (t) Ca " b8 made emp,0yi "9 the decimated and ^slated I and Q signal 
SSUSLKT re ' at,0nSh,p may be expressed as 9 < k > ^ the k-th decimated sample period L The 
phase samples 9(k) are computed from S e (k) via the operation 
S(k) = tan- (Q E (k)/ l E (k)) 
= tan - ' ( sinA«iK/ cosAo»k) 
= tan - ' ( tan(Auk)) 
as = A«k (13) 

where the time interval of k is understood to be the decimated period. 
Three problems arise in calculating tan- (Q E (k) / l E (k) ) in (13)- 

(k) 1 0.' DiViSi ° n ° f ° E(k) ^ ' E(k) 00 3 fraCti ° nal Rxed P ° int maChine is limited to tne case where WO > Qfi- 

lookup 2 ' tan "' (X> iS an inV9rSe transcendental functi °n requiring either a Taylor series expansion or a table- 

, l™ 6 valid ra "9 e of tan "' ( Q e(k) / WO ) extends beyond i180* after phase unwrappina and hence 
the tan- function must be scaled appropriately for a fixed point machine unw 'aPP'ng and hence 

All three problems may be solved by: 

(a) Using a small (256 word) ROM arctangent table over the range (0.45* ]. To use the table a short 
(8 msfruction cycle, division of Q E (k) by l E (k) is required to yield an 8 b* ROM table Ldress Nofe that Zer 
the arctangent ROM table range, the inequality l E (k) > Q E (k) * 0 is preserved 

(b) Extending the range of the arctangent table to £180 ' by 
(i) preserving the signs of Q E (k) and l E (k): 

so (ii) noting whether | l E (k)| > | Q E (k)| or vice-versa; and 

eachoftneVSranf; 600 "^^ "* ° f <° compute the arctangent function in 

ran™ S? !v 6 R ° r M 135,8 Va ' UeS tQ the ,ar 9 est possib,e ra nge of the unwrapped phase e 0 (k) This 
55 7S? ] ^[maximum al.owab.e frequency error (Hz) • duration of 1 time slot (sec) • 360' * 

« 1 GHZ ' - — ™ frequency error is 

t[(2.5kHz)(0.58 msec) 360* + 180"] 
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= :702" = rV. 

All phase values must therefore be scaled by V = 702* to prevent overflow in the phase unwrapping 
algorithm. 

5 

Construction of the phase trajectory from the phase samples. 

Provided that the s.gns of l 6 (k) and Q e (k) are preserved, tan- ( Q E (k) i e (k) ) is defined only over the 
,o .nternal [-». »r]. or equivalently only over the .nterval [-180* . 180 ]. Since absolute phase .s not known in 

the frequency control process, any g.ven sample of S 6 (k) will yield a phase sample d(k> w.th.n th.s range. 

For example, if the output of the phase computation of tf(k) is used directly as the phase traiectory. and if 8- 

(1) = 179 * 3(2) = -179*. 3(3) = -177* then determination of the instantaneous frequency offset via 

time differentiation of the phase trajectory of 3<k) would yield ambiguous results due to the apparent -358 
•5 phase jump between 3(1) and 3(2). In fact, if the phase were "unwrapped" in this example by allowing the 

phase samples to take on values outside the range ^180* . it becomes obvious that the phase samples are 

increasing by +2* every phase sample. 

In order to unwrap the phase samples.. a phase unwrap process 147 is employed in the DSP of tne 

present invention. Conceptually, the unwrapping process for both an increasing phase trajectory and a 
so decreasing phase trajectory is shown in Fig. 2 and .s accomplished by the following steps wh.ch are 

performed by the DSP 1 1 1 of the present invention: 

Initialization: 8(k), k=l,..., N from Step 2 

as Previous = 8(1) ' 

Sumphase = 0 
Threshold = 90* 

DO k = 2 to N 

Current = 8(k) + Sumphase 
IF ( I (Current - Previous)! < Threshold) THEN 
0 o (k) 9(k) 

Previous = 8 0 (k) (No phase jump) 
ELSE IF ((Current-Previous )< - Threshold) THEN 
Sumphase = Sumphase + 360* 
40 8 0 (k) = Current + 360* 

Previous = 8 0 (k) (-360* phase jump) 

ELSE 

Sumphase = Sumphase - 360* 
8 0 (k) = Current - 360* 
Previous = 8 0 (k) (+360* phase jump) 

END IF END DO 

The array of unwrapped phase points 9 0 (k). ke[1.N] are employed in subsequent computation of the 
instantaneous frequency offset. Note that the unwrapped phase samples e c (k) are just the < values A«k given 
by (13) where Awk may now take on any value (including values outside the interval [-180 . 180 ]). 
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Computation of the instantaneous frequency of the frequency-translated signals from the time-derivative of the 

phase trajectory. 
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mJhi h ! , 8b !T!f ,°' 3 7 ,?° iSe ' quanti2ation - frequency instability, the instantaneous frequency error 
may be calculated from (13) using just the difference between any two phase samples In practice al o 

m • T eaSUremen ; °™ f ™ -P'°Ving just two'samp.es. Instead 

linear f t to the phase trajectory ,s employed. The slope of such a least-squares linear fit curve is then 
sZSZTL ° the n in t Sta t ntaneous f ^ uenc V. Prom M. Schwartz and L. Shaw Signal Processing Discrete 
SSl S^S ^r 3 " 0 "; McGr3W - Hi "' ,975 - PP- 14 -' 5 - *• s^sToflhel&sCTaS 
Perform d by function \ 49 ^ ^ e<m " ^ * ^ * ,0 "° Win9 
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N 

X ne 0 (n) 
n=-N 

N 



(14) 



n=-N 
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f - Mi? ha !^r? S ° f d69reeS - ThS instan ^neous frequency error is therefore given by 

ferror(Hz) = dfl/dt = S/360N • L*270.833x10 3 (15) ' 

IfMro mSiS',?"?''" 9 ""f ' emplo ' M "» 0,6 «"> inverters 107 and 109: i.e.. the A/O sampling rat. 

m^sss r src — - — "• « ~ 

• W,th K1 proper se ,ec «0" or the loop gain constant « and knowledge of the VCXO Hzvolt gain constant it 
A summary of these three simulations is given below. 



Simulation 


M 


N 


2N + 1 


Absolute 










Frequency 










Error 


1. 


1 


398 


797 


1Hz 


2. 


20 


19 


39 


11Hz 


3. 


40 


9 


19 


27Hz 
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Since a conventional 8 bit O A converter 155 is employed in the frequency control loop and the stability of 
VCXO 1 15 is :2.5ppm at 1 GHz .mply.ng an error step size of 19 5 Hz. the values of N and M employed in 
simulation 2 w.ll result in a frequency measurement error of approximately 12 LSB o the D.A 155 
referenced to the VCXO 115 stability. The result of simulation 2 is shown in Fig. 3 .n which V c <t) .s plotted 

as a function of frequency offset. 

Note that once N is selected, the denominator in (14) can be pre-computed- For exped.ency. the 
denominator used in <14> can be selected to speed the frequency measurement process by selecting the 



'0 



20 



N 

S' s 2r n 2 S/W 
n=-N 

Us.ng the above-defined constants in (15). the expression for the instantaneous frequency error in terms of 
the measured parameter S is 

f error (Hz) = SY360-N -L • 270.8333 x 103 • WV 
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N 

n=-N 

In the preferred embodiment, 
M = 20 W = 512 

L = 8 V = 2048 

N = 19 N 

Sn 2 = 4940 
n=-N 



4 and thus 

flNSTtH 2 ) 363 ' 875 - 8 ' • 
n=-N 

n=-N 

36 and thus 

ferror (Hz) = 63,875-S' , 



55 A system which would employ a receiver employing the present invention is that which is specified for 
the Pan-European System in the aforementioned GSM documents. A timing diagram shown in Fig. 4 relates 
the burst TDMA transmissions to the frequency control of the present invention. In the exemplary 
transmission shown in Fig. 4 one transmission frame and part of a second transmission frame illustrated as 
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blocks of information transmission conveyed as frequency excursions about a center carrier frequency (U 
Each of these information blocks are conventional TO MA timeslots. each timeslot conveying information to a 
selected different receiver. In the preferred embodiment each frame lasts for a duration of time equal to 
4.6ms. and each timeslot lasts for a time duration of 0.58ms. 
s In order to convey a frequency correction signal from the transmitter to the various receivers the 

preferred embod.ment employs a distinct carrier frequency shift during a predeterm.ned timeslot. such as 
that shown at 401. (This carrier shift has come to be called a "pure sine wave"(PSW)). The carrier shift is 
equal to 13MHz/,92 - 67.708kHz (within a to.erance range of *- O.OSppm = * 0.003»Hz). Regu.ar TDMA 
commun.cat.ons may occur in subsequent timeslots. as indicated. The PSW need not be transmitted during 
w each frame; the preferred embodiment transmits the PSW at a rate of approximately 21 times per second 

The frequency control of the present invention processes the PSW as described previously Upon 
rece.pt of the PSW, the curve relating frequency offset (in Hertz) to the required frequency control signal (in 
m-ILvolts ,s interrogated based on the value of frequency offset calculated by the phase trajectory 
computation. The s.ngle value of control signal <V C (t)) corresponding to the frequency offset is determined 
This singular control signal value is input to the VCXO 115 as a single step that places the VCXO frequency 
within 0.03 parts per million of the carrier frequency of the transmitted TDMA signal. This singular control 
signal value is maintained for the duration of the frames of signal until the next PSW is received In this 
way a rapid and accurate correction of frequency between the transmitted carrier frequency and the 
received frequency is obtained. nuw^y «mu me 
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Claims 



as 



1. A frequency control apparatus for a burst-mode radio communications system employing a fre- 
quency correction signal transmitted as a burst to enable correction of frequency differences between the 
frequency of the rad.o carrier of a subsequent communication burst and the frequency of reception by a 
radio receiver, the frequency control apparatus characterized by- 
means, including a variable frequency oscillator, for receiving a frequency correction signal burst and for 
Processing a communication burst having a frequency difference between said variable frequency oscillator 
30 and the radio carrier frequency of said communication burst; and 

means coupled to said means for receiving said frequency correction signal burst, for generating a singular 
control signal value and for applying said singular control signal value to said variable frequency osciXator 
to correct said frequency difference. 

« f rt H^\ freq ? ^ • Cy , J T tr0, ^P 3 ^ 03 in accordance with claim 1 wherein said means for generating is 
35 further charactered by means for calculating at least one phase magnitude related to said frequency 
difference and wherein said means for processing a communication burst is further characterized by means 
Snafcs" " 9 3 frequency action signal burst into I and Q quadrature frequency correction burst 

.o ^ESSSJg* apparatus in accordance with c,aim 2 wherein said means ,or generatin9 is 

z= L 0 ; at ,east some of said ' and ° quadrature freq — y ^ 

means for sampling said I and Q quadrature frequency correction burst signals to generate I and Q 
correction signal samples; 

means for generating a phase trajectory from at least some of said I and Q correction signal samples- and 
means for compensating a direct current offset. 

4. A frequency control apparatus in accordance with claim 3 wherein said means for calculating at least 
one phase magnitude » further characterized by means for calculating the arctangent of the quotient of a Q 

diV !? ed by a " ' COrreCti ° n Si9na ' Samp,e and wherein s * d ™™ gene at ng a 
phase trajectory ,s further characterized by means for calculating a least squares linear fit cuXe of a first 
.ca.cu.at.on of an arctangent to a second ca.cu.ation of an arctangent, thereby producing a frequency 

?™!hT US "^h Si9 " a ' Va,U6 re ' ati0nShiP fr ° m WhiCh 0ne sinQU,ar control signal value for each 
frequency difference is derived. 

f U rtht'r^h! 1 ? a q ,t enCy H C h 0ntr0, apP f aratus in accordance with claim 1 wherein said means for generating is 
ZZ lZTnT bV m6anS forproducin 9 a f ^ uen <=y difference versus control signa. value relationship 
l° m Wh,ch one s,n 9 u,ar contro ' s, 9 nal val "e for each frequency difference is derived and wherein saiS 

ZnvZ r^T 9 aPP,yin9 S3id S,n9U,ar COntr0 ' Si9nal va,ue to c °™*t »* frequency difference 
applies sa,d Singular control signal value within one frame of a plurality of communication bursts 
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6. An automatic frequency adjustment circuit for a radio characterized by: 
oscillator means for producing a signal at a predetermined nominal frequency. 

means for receiving a communications signal at substantially the frequency of the oscillator means, and a 
reference signal having a fixed frequency offset to the communications signal frequency, and 
5 digital means responsive to the reference signal to produce a correction signal to the oscillator means, to 
cause the oscillator means to produce a signal at the frequency of the communications signal. 

7. An automatic frequency adjustment circuit as set forth in claim 6 wherein the means for producing a 
correction signal is characterized by means for digitally generating a calibrated S-curve of a correction 
value as a function of the frequency offset and which is substantially insensitive to amplitude variations of 

jo the reference signal. 

8. An automatic frequency adjustment circuit as set forth in claim 6 wherein the correction signal 
causes the oscillator means to produce a signal within one part per million of the frequency of the 
communications signal in one step within ten system frames. 

9. A method of frequency control for a receiver in a burst-mode radio communications system 
'5 employing a frequency correction signal transmitted as a burst to enable correction of frequency differences 

between the frequency of the radio carrier of a subsequent communication burst and the frequency of 
reception by the radio receiver, the method of frequency control characterized by the steps of: 
receiving a frequency correction signal burst and processing a communication burst having a frequency 
difference between a variable frequency oscillator and the radio carrier frequency of said communication 
20 burst: and 

generating a singular control signal value in response to said receiving of said frequency correction signal 
burst and for applying said singular control signal value to said variable frequency oscillator to correct said 
frequency difference. 

10. A method in accordance with claim 9 wherein said processing a communication burst step is further 
25 characterized by the step of converting said frequency correction signal burst into I and Q quadrature 

frequency correction burst signals wherein said calculating at least one phase magnitude step is further 
characterized by the steps of calculating the arctangent of the quotient of a Q correction signal sample 
divided by an I correction signal sample, and 

calculating a least squares linear fit curve of a first calculation of an arctangent to a second calculation of an 
30 arctangent,, thereby producing a frequency difference versus control signal value relationship from which 
one singular control signal value for each frequency difference is derived: and 
wherein said generating step is further characterized by the steps of: 

(a> calculating at least one phase magnitude related to said frequency difference. 

(b) translating the frequency of at least some of said I and Q quadrature frequency correction burst 
3S signals by a predetermined frequency, 

(c) sampling said I and Q quadrature frequency correction burst signals to generate I and Q 
correction signal samples. 

(d) generating a phase trajectory from at least some of said I and Q correction signal samples, and 
fe) compensating a direct current offset. 
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